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Abstract. Core-collapse supernova neutrinos undergo a variety of phenomena when they 
travel from the high neutrino density region and large matter densities to the Earth. We per- 
form analytical calculations of the supernova neutrino fluxes including collective effects due 
to the neutrino-neutrino interactions, the Mikheev-Smirnov-Wolfenstein (MSW) effect due to 
the neutrino interactions with the background matter and decoherence of the wave packets 
as they propagate in space. We predict the numbers of one- and two-neutron charged and 
neutral-current electron- neutrino scattering on lead events. We show that, due to the energy 
thresholds, the ratios of one- to two-neutron events are sensitive to the pinching parameters 
of neutrino fluxes at the neutrinosphere, almost independently of the presently unknown 
neutrino properties. Besides, such events have an interesting sensitivity to the spectral split 
features that depend upon the presence/absence of energy equipartition among neutrino 
flavors. Our calculations show that a lead-based observatory like the Helium And Lead Ob- 
servatory (HALO) has the potential to pin down important characteristics of the neutrino 
fluxes at the neutrinosphere, and provide us with information on the neutrino transport in 
the supernova core. 



1 Introduction 



Parallel to the developments occurring in core-collapse supernova simulations [1], the degree 
of complexity in understanding of neutrino flavor conversion in dense matter has enormously 
increased in the last years. It was pointed out some time ago [2, 3] that the neutrino- 
neutrino interaction could trigger new mechanisms for flavor conversion since the neutrino 
number density within the supernova is very large. The ongoing investigations focus on 
unraveling these flavor conversion phenomena, on the underlying mechanisms and/or general 
properties of the associated non-linear Hamiltonian and the phenomenological impact. It is 
now understood that in the inner region, near the neutrinospheres, flavor conversion is first 
frozen by synchronization due to the strength of the neutrino-neutrino interactions [4, 5]. 
Then bipolar oscillation occur [5, 6], ending up in spectral splits and swaps [4, 5, 7-10]. Such 
collective effects can be suppressed by decoherence coming from large matter densities, as 
pointed out in e.g. [11]), and recently investigated using one-dimensional matter density 
profiles from supernova simulations [12]. The presence of invariants of the Hamiltonian 
describing the neutrino evolution has been identified in several works, e.g. [7, 13]. This 
aspect has been further investigated in a recent interesting work [14] where it has been 
shown the presence of invariants not only for the effective one-body mean-field description 
but for the many-body Hamiltonian (when the matter term is not included). The impact 
of the matter phase (and its derivative) - coming from the presence of complex off-diagonal 
contributions from the neutrino-neutrino interaction term - on the neutrino evolution in the 
'matter' basis Hamiltonian is pointed out in [15]. In particular, the start of bipolar oscillations 
is shown to be associated with a rapid growth of the phase derivative. On the other hand 
bipolar oscillations can be understood as the motion of a flavor pendulum [4] and gyroscopic 
pendulum [6, 16]. By building up a self-consistent adiabatic solution in a comoving frame, in 
[7, 17] the authors discuss that the spectral split is a regular precession phenomenon. In [18] 
the spectral split phenomenon is put in relation with a magnetic resonance phenomenon. By 
realizing a full numerical calculation in two-flavors and the single-angle approximation it is 
shown that the neutrino energies and the location where the magnetic resonance condition 
is satisfied are the same as where the spectral split occurs. Putting the spectral split in 
connection with the magnetic resonance phenomenon points out the way to calculate the 
resonance condition in the flavor basis directly, without going to the comoving frame. In [16] 
both precession and nutation driven resonance conditions governing the spectral split are 
discussed. Besides several works have started investigating the phenomenological implications 
of collective effects. The possibility to pin down the neutrino hierarchy and third neutrino 
mixing angle using shadowing from Earth matter effects is studied in [19] or the positron 
time signal associated with inverse beta decay in [20]. The impact on the diffuse supernova 
neutrino background is studied in [21, 22]. While important progress has been made in 
the last years, several aspects still need a full understanding as well as finally assessing the 
phenomenological impact of collective effects in a future core-collapse supernova signal. 

When the neutrino interactions with ordinary matter starts dominating, the Mikheev- 
Smirnov-Wolfenstein (MSW) effect induces efficient flavor conversion if the matter density 
profile is adiabatic [23, 24]. However the arrival of the (front and reverse) shocks in the 
outer layers of the star can leave an imprint on the neutrino signal [25, 26], render the 
conversion non-adiabatic, and induce multiple resonances and phase effects [27, 28]. A full 
calculation of the neutrino amplitudes up to the star surface, using matter density profiles 
from hydrodynamical calculations, has shown that no extra interference effects occur from 
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the combination of the two regions in iron core-collapse supernovae [20]. At the end of their 
trip, decoherence [29] and Earth matter effects [30] impact the neutrino wave packets before 
reaching an observatory on Earth. 

The details of the neutrino spectra at the surface of the star depend upon unknown 
neutrino parameters - in particular the mass hierarchy and the third neutrino mixing angle 
- and upon the "primary" neutrino fluxes (at the neutrinosphere) . In addition, another 
important unknown property is the possible existence of CP violation in the leptonic sector. 
Possible effects coming from leptonic CP violation in media are being explored. First studied 
in [31], it has recently been shown that the Dirac CP violating phase can modify the neutrino 
fluxes in media, both in the supernova environment [32, 33] and in the Early Universe at 
the epoch of Big-Bang nucleosynthesis [34]. Such effects arise for example because of the 
different interaction with matter that muon and tau neutrinos have, due both to radiative 
corrections within the Standard model and to physics beyond the Standard model (e.g. 
Flavor- Changing Neutral currents). The features of the neutrino fluxes at the neutrinosphere 
encode information on various aspects of the supernova dynamics, including the microscopic 
processes determining the neutrino transport within the supernova core and the equation of 
state of the neutron star (for the late time cooling phase). For example, simulations show 
that, in the case of very massive stars leading to failed supernovae and black holes, larger 
average neutrino energies can be achieved compared to the case of supernovae turning to 
neutron stars [35]. 

The neutrino spectra appear to be rather well fitted using Fermi-Dirac or power-law 
distributions [36]. Beside the average energy these spectra are characterized by the second 
moment of the distribution related to the pinching parameter. Supernova simulations suggest 
non-zero values for this parameter [36, 37]. While a hierarchy of average neutrino energies 
is expected due to the difference of the electron neutrino, anti- neutrino and muon and tau 
(anti)neutrino cross sections with matter from the standard model at tree level, the values of 
average energies themselves are subject to variations depending on the supernova simulations 
(see [36] and references therein) . Furthermore the presence of an equipartition of the gravi- 
tational energy emitted during the supernova explosion among the flavors is one of the open 
issues. Indeed, after the prompt burst, the electron-type neutrino luminosities can deviate 
from the non-electron-type (anti)neutrino luminosities within a factor of two, depending on 
the simulations and the phase of the explosion (see e.g. table 7.3 of ref. [38], table 6 in [36] 
and [39]). This also produces qualitatively different neutrino spectra after collective effects 
depending on the primary neutrino flux parameters [40]. 

Several running neutrino detectors are waiting for the next core-collapse supernova 
explosion, while new large size observatories based on various technologies are under study 
[41]. The latter would look for proton decay, CP violation [42] and (extra-)galactic supernova 
neutrinos having, in particular, a discovery reach for the diffuse supernova neutrino back- 
ground [43] . While most of such observatories are sensitive to electron anti- neutrinos through 
scattering on protons, the detection of electron neutrinos exploits scattering on nuclei. 

A new detector is currently under construction at SNOLAB: the Helium And Lead 
Observatory (HALO) [44]. This dedicated supernova neutrino detector is able to observe 
the neutrons emitted from electron-neutrino scattering on lead from charged- and neutral- 
current events [44] . While the idea was first proposed a long time ago [45] different possible 
realizations have been discussed in the literature [46-49]. In [50] it is shown that, among 
different possible nuclei as neutrino targets, measuring one- and two-neutron events on lead 
is particularly attractive to extract information on the supernova neutrino temperatures. A 
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subsequent work [51] has shown that these detection channels are very sensitive to large/small 
values of the still unknown third neutrino mixing angle ^13 and that the possible detection 
of the outgoing electron energy in the detector could help extracting information on the 
pinching parameters of the primary neutrino fluxes. 

The goal of this paper is to investigate the information that can be extracted on the 
characteristics of the neutrino spectra at the neutrinosphere with a detector like HALO, in 
spite of the existing uncertainties coming from the unknown neutrino properties and from the 
different supernova simulations. With this aim we perform new predictions of the charged- 
and neutral-current supernova neutrino signal expected in HALO. The lead-based supernova 
neutrino detector currently built is made of solid lead so that no electron can be measured, 
at variance with some of the detector designs previously considered. Having such a design 
in mind, we perform a detailed analysis of the information that can be extracted from the 
measurement of one- and two-neutron events, if an iron core-collapse supernova explodes in 
our galaxy. Our calculations include for the first time the collective effects due to the neutrino- 
neutrino interactions, the MSW effect in the outer layers of the star and the decoherence 
introduced by the traveling up to the Earth. 

The manuscript is structured as follows. In Section H we present the formalism and our 
analytical treatment of the neutrino propagation, from the neutrinosphere up to a detector 
on the Earth. Section III describes the numerical results on the neutrino fluxes and the 
expected number of events in a lead-based observatory like HALO. The potential to unravel 
information on the neutrino spectra is analyzed. Section IV is a conclusion. 

2 The neutrino evolution from production to detection: the formaHsm 

In this work we will follow the neutrino flavor evolution from the neutrinosphere up to the 
Earth. We evolve the neutrino probabilities. The assumption of propagating the probabilities 
and not the amplitudes is sufficient for our purposes since we do not aim at capturing any 
phase effects, like those arising from the Dirac CP violation phase [32, 33] or from the presence 
of multiple resonances [20, 27, 28]. Possible effects arising from the radiative corrections 
inducing the V^r potential will not be considered here. Note that these might arise in the 
case of very large matter densities [52] or of non-standard contributions from SUSY models 



Therefore, non-electron-type (anti)neutrino fluxes are assumed to be equal and it is 
convenient to work in the I'e, t^x^ basis, which corresponds to mass eigenstate basis at high 
background matter densities close to the neutrino production region [54] . This basis is related 
to the i'e,i^fi, basis through a rotation along the second mixing angle 623- 



such that X- and y-type (anti)neutrino fluxes are equal to other non-electron-type fluxes. 
Following the above assumptions the observable final flavor neutrino fluxes at the Earth 



[53]. 



l^x = COs6'23t'Ai - Sin6'23^'r , 

lyy = sin 6*23 f/x + cos 6*23 i^r , 



(2.1) 



F = (yFy^^Fjj^, Fjy^) can be computed as follows: 




(2.2) 
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where F^^ , F^^ , F^^ are the primary neutrino fluxes at the neutrinosphere, P^u is the neutrino- 
neutrino interaction swapping probabihty matrix, Pmsw is the MSW conversion probabihty 
matrix and A takes into account the decoherence of the neutrino wave packets after neutrinos 
exit the surface. In this notation it is understood that F^,^ = 1 / (Anr'^) F,^^ {£ = e,//, r), where 
r is the distance to the supernova. 

Our eq. (2.2) is assuming factorized dynamics. This hypothesis is based on the fact 
that, within the MSW region, the high and low resonance are well separated in iron core- 
collapse supernova [55] . In [20] a complete propagation of the neutrino amplitudes from the 
neutrinosphere to the star surface has been performed and has shown that it is indeed a 
good approximation to factorize the dynamics between the inner region, where the neutrino- 
neutrino interaction effects are important, and the outer layers, where the MSW resonances 
occur. 

Let us give the explicit expression of each matrix in eq. (2.2). For the collective flavor 
conversion effects in the neutrino-neutrino interaction region, it is sufficient to consider the 
following swapping probability matrix: 



where Pu = P{vi — )• f^) = 1 — Pex — Pey (^ = e, x, y) corresponds to neutrino survival 
probability, while Pex = P{^e ^ ^x) and Pey = P{ve ^ T^y) describe the electron neutrino 
swapping with x- and y-type neutrinos, respectively, such that Pu + Pex + Pey = 1- For 
antineutrinos one can use the same form as above by replacing vi with v^. 

After the collective effects have ceased, the neutrinos enter the MSW region, in which the 
background matter can significantly modify neutrino flavor evolution [23, 24]. One can con- 
struct a matrix consisting of all the possible probabilities that original flavor (anti)neutrino 
state f£ arrives at the surface of the star as a mass eigenstate vi: Pa = P{yi — )• Vi) {£ = 
e,x,y;i = 1,2,3). In normal mass hierarchy (NMH) the MSW level crossing probability 
matrix for neutrinos is given by [55] 



where Ph and Pl are the high and low resonance level-crossing probabilities, respectively. 
The corresponding matrix in inverted mass hierarchy (IMH) is obtained from above by setting 
Ph = 1. In the case of antineutrinos the crossing probability matrix Pmsw in NMH is a unit 
matrix since no resonances occur in the antineutrino sector in this case. The antineutrino 
conversion matrix in IMH is obtained from eq. (2.4) by replacing Pl hy 1 — Pl. Note that 
the high resonance crossing probabilities are equal in the neutrino and antineutrino sectors 



The neutrino flavor evolution depends on how the resonances are crossed. The crossing 
is completely adiabatic {Pl,h = 0) if neutrino remains on its original mass eigenstate and 
completely non-adiabatic {Pl,h = 1) if it jumps to the other energy eigenstate. The adia- 
baticity of the resonances in turn depends on neutrino mixing parameters and the matter 




(2.3) 




PhPl 1-Pl {1-Ph)Pl 
Ph{1-Pl) Pl (1-Ph)(1-Pl) 
1-Ph Ph 



(2.4) 



{Ph = Ph). 
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density profile [55]. Note that the presence of large matter density fluctuations can break the 
H- and L-factorization as discussed in [56, 57]. While in the present work we do not include 
either shock wave or turbulence effects a discussion will be made of how our results might 
be changed if one implements supernova matter density profiles with a (front and reverse) 
shock [20, 27, 28] and matter density fluctuations [57, 58]. 

Once the (anti) neutrinos have reached the surface of the star as mass eigenstates'^] (i = 
1,2,3), they travel up to the Earth and are detected as flavor eigenstates vl {I = e,/U,r). 
For typical distances between a supernova in our Galaxy and the Earth the mass eigenstate 
wave packets decohere [55] so that no extra effect coming from the oscillations occur. This 
is taken into account by the use of the matrix: 

A = I \U^2? \U^^? , (2.5) 

\Url? \Ur2? \Ur3\y 

consisting of all the probabilities of a certain flavor neutrino being in a given mass eigenstate 
i^i- \Uii\^ = {i = e, jjL,T\i = 1,2,3). Ua are the elements of the Pontecorvo-Maki- 

Nakagawa-Sakata (PMNS) mixing matrix C/pMNS which is the unitary rotation matrix relating 
the flavor and the mass basis. Considering three neutrino flavors, such a matrix depends on 
three mixing angles 0i2, ^13, ^23; one Dirac and two (possible) Majorana phases. The standard 
parametrization of this mixing matrix can be found e.g. in [59]. 

Finally, let us discuss the primary neutrino fluxes at the neutrinospheres, in eq. (2.2), 
whose features reflect the underlying neutrino microscopic processes in the neutrino transport. 
According to current supernova simulations, that include neutrino transport by solving the 
Boltzmann equation for particles without mixings (see e.g. [60]), the neutrino fluxes can be 
rather well parametrized either by Fermi-Dirac or modified power-law distributions. As a 
consequence, can be defined as 

F'AEu) = ^<PiE.) , (2.6) 

where L,y is the neutrino luminosity, {E^) the average neutrino energy and <j){Ey) the neutrino 
energy spectrum at the neutrinosphere. 

For normalized (J dEcl){E) = l) "pinched" Fermi-Dirac distribution the energy spec- 
trum can be written as 

1 E^ 

T^F2{riy) exp {E^l% - r?^) + 1 

where is the neutrino temperature, r]y is the degeneracy or pinching parameter and the 
average neutrino energy can be expressed as 

(E»=T„||). (2.8, 

The functions F2{r]u) and F2,{r]^) are the Fermi integrals: 

FuM = / , . (2.9) 

7o exp (x -r]u) + l 

The relation between r]y and T^j is plotted in 2 using different values for {Eu). 



-6- 



In [36] it is argued that the power-law parametrization is more flexible at describing the 
high-energy tail of the neutrino spectra in comparison to the above Fermi-Dirac distribution. 
For normalized modified power law or alpha-fit distribution the energy spectrum is given by 



where au characterizes the pinching and 



(2.10) 



{E,)r{a, + 1) 

is the normalization constant, in which T{p) is the Euler Gamma function.^ 

In the future, it would be very valuable to extract, from supernova simulations, more 
information on the pinching parameter (and on a possible non-thermal component) of the 
neutrino energy distributions. It is indeed one of the goals of the present work to emphasize 
the importance of the tail in determining the response of detectors having neutrino detection 
channels with high energy thresholds such as the two-neutron one in HALO or neutrino- 
carbon scattering in a scintillator detector like LENA [61]. Obviously, if one sticks to the 
above-mentioned parametrizations (2.6-2.10), as we do in the following, the primary super- 
nova neutrino fiuxes are fully determined if the neutrino luminosity, the average energy and 
the pinching parameter are identified. 



3 Numerical results 

The goal here is to present predictions for the charged-current and neutral current events 
in a lead based detector like the Helium And Lead Observatory (HALO). The detector is 
currently under construction and should start to be operating soon [62]. In phase-I HALO 
is made of 79 tons of lead while in phase-II the mass is planned to be increased to 1 kton. 
Such a detector exploits the following detection channels: 

Ue + ^°^Pb ^ 2°^Bi + 2n + e" 

+ 208pb ^ 206p^j + 2n + 

where the neutrons are detected using ^He counters as done for the SNO experiment [63]. 
HALO will measure neutral- and charged-current events without the capability of distinguish- 
ing them since the outgoing electrons are not detected. In phase-I the one-neutron (In-) and 
two-neutron (2n-) detection efficiencies are about 50 % and 25 % (for seeing both neutrons^), 
respectively, while the detector will have a good time resolution of about 30 ms [64]. In the 
case of a supernova explosion, HALO will have the capability of following the different phases 
of an explosion, from accretion phase to the cooling of the proto-neutron star. The specificity 
of measuring In- and 2n-events gives to such a detector an interesting physics potential, as 
we will discuss in the following. 

'^We call the parameters rji, and simply as pinching parameters since their variation directly either 
suppresses or enhances the high energy tails of the distributions. 

^If both neutrons are captured efficiency is nearly f 00 % due to their small separation in time. Additionally, 
in practice, some 2n-events can give f n-signal if the other neutron escapes the detector. Moreover, secondary 
neutron production can be seen as a 2n-signal, although unlikely. Such corrections should be determined by 
measuring the detector response using HALO phase-I. 
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Neutrino temperature, Ty [MeV] 

Figure 1. Neutrino degeneracy parameter 77^ as a function of neutrino temperature with 
different average energies (see eq. (2.8)). Neutrino average energies are from left to right: 
10,13,15,18,21,25 MeV. 



3.1 Choice of the parameters 

We consider that a core-collapse supernova explosion occurs in our galaxy, at about 10 kpc 
from the Earth. For the neutrino fluxes at the neutrinospheres in eq. (2.2), both Fermi- 
Dirac eqs. (2.7)-(2.9) and a power-law spectra eq. (2.10) have been investigated. As an 
example, we consider luminosity values corresponding to later cooling phase of the proto- 
neutron star. In order to test whether the energy equipartition is satisfied or not during 
this time period, we consider here two possibilities: an equal luminosity case with L{_) = 

1/6 X 10^^ erg s~^ and an extreme case of Lj^^ = 2Lj^^ with L^^ = Lp^ = 1 x 10^^ erg s~^, such 
that the total neutrino luminosity LJ^""" = L^^ + Lp^ + 4:L^^ = 10^^ erg s^^ is the same in 
both cases in order to be able to compare them.'^ Here(after) we refer to the all non-electron 
type neutrinos simply by Vx (unless otherwise mentioned). 

The primary neutrino average energies satisfy the hierarchy {E^J < {Ep^) < {Ep^), 
as expected from the difference among the electron (anti) neutrino and the muon and tau 
(anti)neutrino interaction cross sections with matter, and also from the fact that the medium 
is neutron-rich. In our calculations we have fixed {E^J = 10 MeV, {E^J = 13 MeV and for 
the non-electron-type neutrinos we consider the values {E^_^) = 13,15,18,21 and 25 MeV. 

^As suggested by figure 1 in ref. [40], when 0.5 < L^^/L^^ < 1 there is no quahtative difference in neutrino 
spectra in collective region within the considered range of neutrino average energies. Therefore, the other 
extreme case of Li,^ = O.SLi^^ is qualitatively the same as our equal luminosity case. Although, quantitatively 
our results in numerical section do depend on the exact values of luminosities. 
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Neutrino Energy [MeV] Neutrino Energy [MeV] 

Figure 2. Neutrino energy spectrum at the neutrinosphere. Top left: Power-law (PL) spectrum, 
eq. (2.10), with a^, fixed to value 7 and average energies, {E^), vary from 10 to 25 MeV. Top right: 
PL spectrum with values varying from 2 to 7 for two different average energies (10 and 25 MeV). 
Lower figures: Comparison of the initial neutrino Fermi-Dirac (FD), eq. (2.7) (solid), and power-law 
(PL, dashed lines) energy spectra for two different average energies (10 and 25 MeV). Bottom Left: 
?7i/ for FD and = 2 for PL. Bottom right: rj^ ~ ^ and a^, = 7. 

Indeed although fits to the SN1987A data appear to show rather low energy components 
of the fluxes, we wish to explore the sensitivity of our findings to the variation of such 
input parameters, without sticking to any particular model, having in mind that, e.g. stiffer 
equations of state for the neutron star might give higher average energies [35]. 

In all of the cases considered we have fixed the electron- type parameters to oty^^p^ = 3 
and rjy^^y^ = 1 and vary the non-electron-type parameters model independently: rjy^ G [0,5] 
and ajy^ G [2,7]. Note that in ref. [36] all the neutrino fluxes are found to be characterized 
by pinching parameters with values G [2,5] for power-law, eq. (2.10), or rjy G [0,4] for 
Fermi-Dirac distribution, eq. (2.7). We have tested that varying the electron-type pinching 
parameter values has only a very small or no effect in all of the cases we consider. Currently 
still little information on such parameters can be found in the literature from supernova 
simulations. However it is clear that to properly define neutrino fluxes one needs the first as 
well as the second moment of the distribution [36, 37]. Figure 2 presents the neutrino fluxes 
at the neutrinosphere for the two parametrizations considered. Since the results we obtain 
are qualitatively and quantitatively very similar in the following we show results obtained 
for the power-law spectra only. 

As far as the neutrino masses and mixing are concerned, most of such parameters are now 



-9- 



known, in particular Amj^ = (7.59 ±0.20) x lO'^ eV^ {Aml^l = (2.43 ±0.13) x lO'^ eV^(Ri 
|Am|i|), sin^ 26*12 = 0.87 ± 0.03, sin^ 26*23 > 0.92 and sin^ 26*13 < 0.15 [59]. For the mixing 
angles we have used the values sin^ 2^12 = 0.86, sin^ 2^23 = 0.99 while for the third angle our 
results are obtained either considering that ^13 is very small (0.001) or close to the present 
Chooz limit (0.1). We shall consider both normal (inverse) hierarchy coming from a positive 
(negative) sign for |Am23|. 

Note that the measured values of the solar parameters make the flavor transition at the 
low resonance always adiabatic. That is, in eq. (2.4) Pl = for typical density profiles from 
supernova simulations, while the adiabaticity of the high resonance crossing depends on the 
still unknown third neutrino mixing angle ^13. Our choices for the ^13 imply that all the 
processes are assumed to be completely adiabatic Ph = (large) or non-adiabatic Ph = 1 
(very small ^13). We set the Dirac CP phase to zero. For a discussion of its impact see 
[32, 33]. 

3.2 Computation of neutrino fluxes 

For the P^u matrix eq. (2.3) we parametrize the neutrino fluxes after the region where the 
neutrino-neutrino interactions are important following the results of the full three flavor 
numerical calculations with/without equipartition for the neutrino luminosity [10, 40, 65]. 
These references show that after the collective effects there can be no splits, one low (or 
high) energy split or both low and high energy splits in (anti)neutrino spectra, depending on 
the neutrino flux parameters and mass hierarchy.^ In our considerations this implies that in 
eq. (2.3) we take 

Pe <^<^^ (3.2a) 



0, otherwise 

1, if ^ > Ef^ 
0, otherwise , 



(3.2b) 



where, for simplicity, we have assumed all the splits in energy to be sharp i.e. the swapping 
probabilities are step functions. Ef and Ef^ are the low and high split energies, respectively. 

The above parametrization can be easily modified to correspond other possible cases. 
For instance, it is to be noticed that the recent multi-angle simulations [65] show differences 
in neutrino spectra after the collective effects compared to the single-angle simulations given 
in refs. [10, 40]. General feature of multi-angle effects is to make splits less sharp in energy. In 
IMH ref. [65] shows that considering neutrino spectra with single split (our equal luminosity 
case) multi-angle effects are marginal, but for multiple splits (our Lj^^ = 2Ljy^ example case) 
multi-angle effects can be relevant. Namely in IMH with multiple splits one finds strong 
suppression of (e — x)-swapping. In our case for IMH with L^^ = 2L^^ this multi-angle 
result is obtained by setting in eq. (3.2b) Pex = at all energies, that is, there exists only 
(e — |/)-swap at intermediate energies and at low and high energies spectra remain equal to 
the primary. 

Using eqs. (3.2) and (2.3) one can compute the fluxes after the collective region F' = 
(F' F' F' )'^- 

F' = P^^F^ . (3.3) 



"^The current status for the appearance of the sphts in neutrino energy spectrum and their dependence on 
neutrino luminosities and mass hierarchy can be found e.g. in figure 1 of ref. [40]. 
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Table 1. Neutrino fluxes after the collective region {F^^, F^^, F^^) given by eq. (3.3), depending on 
their energy. The fluxes are given in terms of the neutrino fluxes at the neutrinospheres {F^ ,F^ ,F^ ), 
eq. (2.6), in equal and non-equal luminosity case, both in normal (NMH) and inverted mass hierarchy 
(IMH). Ef and Ef^ are the values of the low and high split energies, respectively (see text). The | 
sign means that the result is the same as above. 

The fluxes are listed in table 1 in the cases of equal and non-equal {L,^^ = 2Ly^) luminosity 
both in normal and inverted mass hierarchy. According to [10, 40], in the computation of 
event numbers we have used the numerical values E^ = 8 MeV and Ef^ = 23 MeV, when 
appropriate (see table 1).^ 

The exact value of the split energy depends on the neutrino mixing and the primary 
neutrino flux parameter values. Considering only two neutrino flavors the equation of motion 
for the system implies the "conservation of flavor lepton number" (see e.g. ref. [66]): 

N,^ - iVp, = D{r) j dE, {El - fI) = const , (3.4) 

where Ny^[Ni)^) is the effective number density of fe(i/e), D[r) includes (flavor independent) 
geometrical dependence and F^^ is given in eq. (2.6). 

In core-collapse supernovae initial z/g number density is larger than that of Ve- Therefore, 
if all the U^s convert to non-electron antineutrinos the conservation of flavor lepton number 
implies that not all of the v^s can convert to non-electron neutrinos. This condition is met 
e.g. in the case of equipartition of energy in IMH. Assuming complete Ve -H- Vy swap in the 
antineutrino sector, conservation of flavor lepton number in eq. (3.4) indicates a single split 
energy E'^ for neutrinos: 

/ dE^iF^^ - F^J - (Fl - El)] = dK{Fl-Fl)- dE,.iF^^ 

Jo Jo JE" 

POO 

+ / dE^Fl-El), 
Jo 

where l.h.s. corresponds to initial and r.h.s. to final states, giving 

/ dF.(FO - FOJ = / dF.(FO -El). (3.6) 
Je'' Jo 

We have tested that using eq. (3.6) the value of E^ varies from 6 to 10 MeV when 
changing the primary neutrino flux parameter values (r/ or a^, {E^)) within the ranges given 

^For antineutrinos one can use the same form as above by replacing by f'^ = (F^^, F^_^, F^^)"^ and P^i, 

with P^^ such that Pu.ex,ey = Pu,ex,ey and Ei f^ = E^l^h- 



_pO ) 

^ Vxl 

(3.5) 
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in the previous subsection. Such a modification of has no consequence on the results 
presented because of the location of the In energy threshold. Therefore, we have kept the 
split energy values fixed in all the cases. On the other hand, we have checked that considering 
a similar variation of ±2 MeV for the higher energy split has only a small effect on the numbers 
of events in a Pb detector that we will show. 

The (primed) fluxes after the collective region in eq. (3.3) enter the MSW region. The 
mass eigenstate fluxes F^. which exit the surface, can also be calculated using eqs. (2.2) 
and (2.4): 











= Pmsw 











(3.7) 



As pointed out earlier, the fate of the neutrino fluxes in the MSW region depends on the 
neutrino mass hierarchy and the adiabaticity of the resonance crossings. For neutrinos in 
normal mass hierarchy for small ^13 {Ph = 1) and large ^13 {Ph = 0) the neutrino mass 
fluxes which exit the star are given by 

for small ^13 . 

(3.8) 

for large ^13 . 

In IMH Ph = 1 regardless of the value of ^13. Therefore, in the absence of collective effects 
the case of inverted is degenerate with normal mass hierarchy and small ^13 . 

The observable neutrino fluxes at the Earth are now obtained from the mass eigenstate 
fluxes by using the matrix A in eq. (2.5). That is, for instance, the final electron neutrino 
flux at the Earth is given by 

3 

F,^ = Y,\Uei? ^ (3.9) 

i=l 

where using 6*13 = 0.001(0.1) 

\Uei\^ ~ 0.69(0.68) , 

|C/e2p ^ 0.31(0.31) , (3.10) 

|C/e3|' « 10-^(0.01) , 

giving the weights of the mass eigenstate fluxes which exit the surface. Table 2 shows the 
energy dependence of the mass eigenstate fluxes in terms of the primary neutrino fluxes in all 
the scenarios considered. In figure 3 the corresponding final electron neutrino fluxes (using 
the power law distribution in eq. (2.10) for the primary fluxes) are plotted as a function of 
energy, also showing the dependence of the spectral shape on the value of a^^. The figure 
also presents the neutron emission charged-current Ve — Pb cross sections, to show which part 
of the Ve spectrum detected on Earth is relevant from the point of view of neutrino detection 
on a lead observatory. 
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Table 2. Neutrino mass eigenstate fluxes {F^^ , F^^ , F^^) , eq. (3.7), depending on neutrino energy, 
exiting the star after the MSW resonance region in equal luminosity case. The fluxes are given in 
terms of the neutrino fluxes at the neutrinospheres {F^^, F^^, F^^), eq. (2.6), both in normal (NMH) 
and inverted mass hierarchy (IMH). Notice that in normal mass hierarchy these results depend on 
the value of the third neutrino mixing angle ^13. i?f and are the values of the low and high split 
energies, respectively (see text). The || sign means that the result is the same as above. 




Figure 3. (Color online) Electron neutrino fluxes at Earth, eq. (3.9), (solid lines) as a function 
of energy including v — v interactions, the MSW effect and decoherence. The primary v^^ (blue thin 
dashed) and 1/3; (black and red thin dashed) fluxes at the neutrinosphere, eqs. (2.6) and (2.10), are also 
shown. Black (red) lines correspond to pinching parameter a,j^ = 2 (7). For the primary v,, flux 
= 3. In this figure the primary average energies are fixed as (i?°^) = 10 MeV, — 18 MeV. Upper 

row: equal luminosities, lower row: L,^^ = 2i,y_,. Left panel: inverted mass hierarchy (IMH), middle 
panel: normal mass hierarchy (NMH) with large 6'i3, right panel: NMH with small ^la. Additionally, 
the charge current i^e — Pb one-(thick, solid, orange) and two-neutron (thick, dash-dotted, orange) 
emission cross sections (from ref. [51]) are also shown. 



3.3 Computation of event rates 

Since the detector we will consider does not identify the outgoing lepton, to compute the 
total event rate N"^^"^ we have to add charged-current (CC) and neutral-current (NC) event 
rates 



(3.11) 
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Table 3. Neutrino mass eigenstate fluxes [F^^, F^^, F^^), eq. (3.7), depending on neutrino energy, 
exiting the star after the MSW resonance region but for the L^^ = 2Li,^ case. The fluxes are given in 
terms of the neutrino fluxes at the neutrinospheres (F°^ , , F°^), eq. (2.6), both in normal (NMH) 
and inverted mass hierarchy (IMH). Notice that in normal mass hierarchy these results depend on 
the value of the third neutrino mixing angle ^13. Ef and E\ are the values of the low and high split 
energies, respectively (see text). The || sign means that the result is the same as above. 



where iVi„(2n) refers to ln(2n) event rates. The CC and NC event rates are given by 



iV- 



CC 
ln(2n) 



iVpb / dF^^F,,(F) 



^^ln(2n) 



Pb 



dF 
dF 



(3.12) 



da 



NC,u 
ln(2n) 

dF 



da 



NC,v 
ln{2n) 

dF 



(3.13) 



with A^pb the number of target particles in lead detector^, d(Ti„(2n)/dF the differential cross 
sections and F(_) the final (anti) neutrino flux given by eq. (2.2)7 

As an example, we consider the case of inverted mass hierarchy (single-angle) with 
Ly^ = 2Ly^. According to eqs. (3.9) and (3.12) (results listed in table 3) the charged-current 
event rates are explicitly calculated in terms of primary neutrino fluxes given in eq. (2.6) as 
follows: 



ivrCC 
^^ln(2n) 



pa A CC 

^^ln(2n) 



dF- 



dF 



2 rO 



|2 T?0 



(2n) 



^E ,^ 
dF 

^^ "^ln(2«) 

Es dF 



2 zpO 



|2 



(3.14) 



Uel\'Fl + \U,-2\'Fl + \Ue2\'F[ 



2 770 



|2 77O 



where the weights \Ueif' are given in eq. (3.10). Note that, according to our assumptions, 
FO = FO 

Vy Vx' 

^Here we consider 1 kton lead made of ^°*Pb so that A''pb = 2.9 x lO"^*^. Our results can be easily scaled to 
natural composition. 

'^Notice that, X]f=e u t ^(-) ~ X^f=e u t -^''(-)i that is the NC events are not affected by the flavor trans- 
formations. 
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Figure 4. j/g — Pb charged current total (0 + ln + 2n) cross sections as a function of energy; RPA [51] 
(solid), CRPA [67] (dashed) and RRPA [68] (dash-dotted). 



The neutrino- lead total (0 + In + 2n), In- and 2n-emission cross sections we use are 
from a microscopic calculation based on the Random Phase Approximation (RPA) (see table 
I of ref.[51]).^ These calculations include nuclear matrix elements up to multipolarity 4 and 
assume neutrons emitted with a 100 % probability once the neutrino energy is above the 
In- and the 2n- emission thresholds^. A refinement of such a calculation would require the 
inclusion of all possible decay channels with a statistical code, as done in [67]. 

In order to show the current uncertainty on the neutrino-lead cross sections figure 4 
presents a comparison of the total (0 + In + 2n) neutron emission cross sections with those of 
refs. [67, 68] based on Continuum RPA (CRPA) and Relativistic RPA (RRPA), respectively. 
Figure 5 shows ratios of flux-averaged neutrino-lead cross sections with the neutrino fluxes 
given by the power-law distribution with varying average energies and pinching parameters. 
One can see that, even for low average energies, the uncertainty coming from variances of the 
current theoretical predictions ranges from about 1.2 to 1.4. In the following we will show 
that, still interesting information can be extracted with a lead-based detector, in spite of the 
existent cross section uncertainties. 

In all of our calculations we assume a 100 % detection efficiency, so that our results can 
be scaled easily, depending on the actual neutron detection efficiencies^*^. The results shown 
in the next section correspond to HALO phase-II (1 kton), while those for HALO phase-I 
can be obtained by multiplying the rates by a factor of 0.079.^^ 



The electron anti-neutrino detection channel is strongly suppressed due to Pauli blocking. 
®In [51] it is assumed that the probability for emitting one (two) neutron(s) is 100 % if the neutrino energy 
is above 10 MeV (20 MeV). 

^"^Strictly speaking, detection efficiencies do have some dependence on emitted neutron energy, although it 
is likely to be small. A refinement of the present calculations would require to determine the neutron energy 
spectrum as done in [67] 

^^Note that, if the absolute value of luminosities, detection efficiency, size of the detector or distance to 
supernova are different from what considered here, one can just multiply our results with a suitable factor as 
long as the change does not pose qualitative differences. 
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Figure 5. Contour plots of ratios of flux averaged total (0 + In + 2n) cross sections. Left: RPA [51] 
to CRPA [67]. Contours correspond to ratios from 1.1 to 1.26 (outside to inside) with step of 0.02. 
Right: RPA to RRPA [68], similarly from 1.4 to 1.24 with step of 0.02. 



3.4 Expected events in HALO 

We present new predictions for the In- and 2n-events going beyond the previous calculations 
[50, 51] and include collective effects, assuming an explosion in our galaxy. In particular, the 
results shown are : i) for event rates ; ii) for total number of events ; iii) for ratios of events. 
Since HALO has a good time resolution, the event rates are given in order to make time 
dependent predictions that should be multiplied by the time duration of the corresponding 
phase in the explosion. In such a period the neutrino flux parameters are assumed to be 
fixed. Then we show, as an example, the expected numbers of events during the accretion 
and cooling phases and from the whole explosion. Finally, ratios of one- and two-neutron 
events are given since these turn out to be particularly sensitive to the primary neutrino 
fluxes. Besides they have the advantage that they are free from common normalization 
parameters, like the total luminosity. 

Our goal here is two- fold. First, we want to give updated predictions for the HALO 
project under construction at SNOLAB. Second, we investigate which information on the 
neutrino fluxes at the neutrinosphere and/or on the neutrino properties can be extracted 
with such an observatory. For example, it is a priori not clear that information can be 
extracted from the primary pinching parameter, after the inclusion of the collective effects 
and with the current design of the HALO detector. Indeed it was first pointed out in ref. [51] 
that the measurement of the energy of the outgoing lepton would help constraining the 
primary pinching parameter. As we will discuss, by considering combinations of In- and 
2n-events some information on the primary neutrino fluxes can indeed be obtained. 

Before presenting our results we emphasize that, in order to extract as much information 
as possible from a future measurement, it will be indeed crucial to combine results from the 
HALO detector with those from detectors based on other technologies, like water Cherenkov 
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(e.g. Super-Kamiokande or MEMPHYS), scintillators (Borexino, LVD, KamLAND, LENA) 
or liquid argon (GLACIER) [41]. These detectors will be sensitive to electron anti-neutrinos 
through inverse beta decay and electron neutrinos using scattering on nuclei (oxygen, car- 
bon, argon). In particular, the measurement of neutral current events through scattering on 
electrons or protons [69] can be used to determine the total neutrino luminosity and possi- 
bly also extract the spectra of non-electron type neutrinos. Note also that the number of 
unknowns considered in the present work might be narrowed down in the near future if the 
third neutrino mixing angle is determined by the reactor Double-Chooz, Daya-Bay, RENO 
or T2K and NOi^A experiments [70]. 

Results for In- and 2n-event rates, as a function of non-electron-type primary neutrino 
pinching parameter a,y^ , are given in figure 6 for equal luminosities and in figure 7 for Li,^ = 
2Ly^ , with different non-electron- type primary neutrino average energies {E^^ ) . These results 
show how the variation of the parameters describing the primary non-electron type neutrino 
energy spectrum affects the observable event rates. The results for the normal mass hierarchy 
with large ^13 are degenerate with inverted mass hierarchy. This is because the low energy 
split present in the former case is not visible due to the high neutron emission thresholds in 
lead. It is for the same reason that the prompt Uf> burst occurring at the very early stages of 
a possible future core-collapse supernova explosion will not be observable in such a detector 
if the primary electron neutrino average energies are as low as 10 MeV. Note that, in the 
absence of collective effects, the results of IMH would be the same as for the NMH with small 
^13 shown here. It is clear that the In- and 2n-events alone will not allow the identification 
of all the degenerate physical scenarios depending on either the primary flux parameters 
{{E^^), Oi,^, neutrino luminosities), or neutrino properties (mass hierarchy and the value of 
^13). The neutrino luminosity, that gives the normalization factor of the events shown, can 
be obtained from the measurement of neutral current events in other detectors, as discussed 
above. From figures 6 and 7 it can be seen that In- and 2n-events provide complementary 
information: the In-events are more sensitive to pinching with smaller values of {E^J)^ while 
the 2n-ones are sensitive to high average energy values. 

As far as number of events are concerned, the range of events expected during the 
accretion phase are given in table 4 in different physical scenarios. For this case a time 
duration of 300 ms is assumed and a total luminosity of 18 x 10^^ erg s~^.^^ Therefore, the 
numbers are obtained by multiplying the events given in figure 6 by a factor 5.4 s. It is to 
be noticed that the collective effects can be suppressed by the background matter during 
the accretion phase [12]. Even if this finding is consolidated, our results remain useful since 
IMH would be then degenerate with NMH and small ^13 and in NMH collective effects are 
not present in the neutrino channel. Table 5 presents the total number of events for the 
cooling phase using results in figures 6 and 7 multiplied by a duration time of 10 s. In table 6 
the total numbers of events for the whole explosion are given assuming equal luminosities 
throughout the neutrino emission with the total time-integrated neutrino luminosity fixed at 
the typical value of 3 x 10^'^ erg^''. 

Since In- and 2n-events can be well identified in HALO, it is attractive to consider the 
corresponding ratio (figure 8), which is independent of common normalization factors (such 

^^According to e.g. [39] for 18 solar mass iron-core progenitor the duration of accretion phase is about 
350 ms and the electron-type neutrino luminosities on average about 0.4 x 10^"' erg s~^ (about half of this for 
non-electron-type) . 

^^These numbers are obtained by multiplying the results in figure 6 by a factor of 30 s to obtain the given 
total time integrated luminosity. 
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Figure 6. One- (upper) and two-neutron (lower row) total (charged -I- neutral current) event rates 
as a function of primary non-electron neutrino pinching parameter ai^^ , including v — v interactions, 
the MSW effect and decohercnce. Here is assumed 100 % detection efficiency and equal luminosities, 
erg s ^, in IMH (left) and NMH with small ^13 (right panel). Distance to the supernova 
is taken to be 10 kpc and target mass 1 kton of ^''^Pb. Different lines correspond to different non- 
electron-type primary neutrino average energies (from bottom to top: 13, 15, 18, 21 and 25 MeV). 



as a distance to the supernova and a size of the detector). From figure 8 it is evident that 
these ratios are sensitive to the pinching regardless of the physical scenario i.e. luminosities, 
neutrino mass hierarchy and the value of ^13.^^ Since in HALO, if a supernova explodes, 
one measures a given ratio, figure 8 shows that the measured ratio would allow to identify 
different degenerate combinations of non-electron-type primary neutrino average energies 
and pinching parameters (even without a precise knowledge of the above mentioned physical 
scenarios and common flux parameters). By knowing (or assuming) the non-electron-type 
primary neutrino average energy^^, it would be possible to give tight constraints on pinching. 



^*Notice the curves in the case of equal luminosities are within the ones of IMH (solid) and NMH with small 
013 (dashed) in = 21/^^ case. 

^®In the cases we have been considering the final electron neutrino fluxes are always close to the non- 
electron-type neutrino primary fluxes (see figure 3). More precisely, in the equal luminosity case, both in 
inverted and normal mass hierarchy with large ^13, the final electron neutrino energy spectrum is the same as 
the primary energy spectrum of non-electron type neutrinos due to the complete flavor conversion (in IMH 
after the low energy split). In normal mass hierarchy with small ^13 there are less electron neutrinos at higher 
energies (relevant from the detection point of view) with respect to the primary non-electron type neutrinos. 
In the = 2L^^ case the average energy of the observable electron neutrinos is somewhat smaller than that 
of the primary non-electron type neutrinos in all of the studied cases. Notice that the figure 3 corresponds 
only to the case with {E'l_^ ) — 18 MeV while the exact relation between the observable and the primary Vx 
spectrum obviously depends on the exact value of (-B°^). 
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Figure 7. One- (upper) and two-neutron (lower row) total (charged + neutral current) event rates 
as a function of primary non-electron neutrino pinching parameter a^^ , including v — v interactions, 
the MSW effect and decoherence. The case is of L^^ = IL^^ (Lp^ — L^^ = lO'^^ erg s^^) in IMH 
(left), NMH with large ^13 (middle) and NMH with small 6*13 (right). Different lines correspond to 
different non-electron-type primary neutrino average energies (from bottom to top: 13, 15, 18, 21 and 
25 MeV). We assume 100 % detection efficiency, L™''" = 10^^ erg s~^, distance to the supernova 10 
kpc and target mass 1 kton of ^"^Pb). 
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~ 16 


124 61 
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Table 4. Accretion phase: Total numbers of one- and two-neutron events, Ni^ and N2n, respectively, 
and total number of neutrons emitted (iVi„-|-2A^2n): depending on the value of the pinching parameter 
a^^ and neutrino properties, including v — v interactions, the MSW effect and decoherence. Here is 
assumed equal neutrino luminosities and 100 % detection efficiency. The total (sum of all flavors) 
time- integrated luminosity is taken to be 5.4 x 10^^ erg (see text for details), distance to the supernova 
10 kpc and target mass 1 kton of ^"^Pb. The column with 13 (25) MeV represents the scenario with 
smallest (largest) amount of events. In addition, the range of possible event numbers are shown for 
the intermediate case with 18 MeV. Notice that the results for NMH and large 6*13 are degenerate 
with ones in IMH shown here. 

The more favorable case is the one with equal luminosities since the degenerate combination 
of average energies and pinching parameters would have only a slight dependence on the 
unknown neutrino properties (figure 8). Note that, while in very few cases the event number 
is small, in the majority of cases statistics is sufficiently large so that the statistical error is 
not a limiting factor. This is shown in figure 9(b) where the area above each curve gives the 
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Table 5. Cooling phase: Total numbers of one- and two-neutron events, A^i„ and N2n, respec- 
tively, and total number of neutrons emitted, depending on the pinching parameter a^^^ and neutrino 
properties, including 1/ — v interactions, the MSW effect and decoherence. Here is assumed that 
L^^ = 2L^^ (equal luminosities in parenthesis) with the total time integrated luminosity 10^'^ erg and 
{E^J = 18 [MeV]. The calculations are obtained using 100 % detection efficiency, a distance to the 
supernova 10 kpc and a target mass 1 kton of ^"^Pb. The column with 13 (25) MeV represents the 
scenario with smallest (largest) amount of events. 
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Table 6. Total numbers of events during the explosion (assuming 100 % detection efficiency, distance 
to the supernova 10 kpc and target mass 1 kton of ^"*Pb). As in table 4 but assuming equal neutrino 
luminosities throughout the whole neutrino emission and the total time integrated luminosity 3 x 

relative statistical errors associated with the event ratios of figure 8. 

In figure 10 we give ratios of 2n-events to determine which precision should be achieved 
in measuring the numbers of events to be able to distinguish different physical scenarios. To 
quantify this, First, we assume the neutrino mass hierarchy to be known (inverted or normal 
with small ^13) and present the ratios of the 2n-events for the case with equal luminosities 
over the one with L,^^ = 2L^^. Second, we fix the luminosity (either equal or L^^ = 2L^^) and 
show the ratios of the 2n-events with IMH over NMH with small ^13. While these ratios are 
obviously not observable, we employ them to give an idea of the required precision to be able 
to identify either different neutrino properties or supernova emission cases. From figure 10 
one can clearly see that the difference between the results in different scenarios range from 30 
% up to a factor of 2^^. The largest difference is due to the presence of the high energy split 
in the neutrino fluxes. Note that a detector capable of measuring the neutrino spectra like 
a scintillator or a liquid argon supernova observatory can actually observe the high energy 
split that would further restrict the range of possible primary neutrino flux parameters. 

^®Note that, while we use a difference between L^^ and L,y^ as large as a factor of 2, similar results are 
obtained if such a difference is as low as 30 %. 
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Figure 8. Ratios of one- (-/Vi„) to two-neutron (-/V2„) event rates as a function of primary non- 
electron neutrino pinching parameter a^^ with different primary non-electron-type neutrino average 
energies. Predictions include — interactions, the MSW effect and decoherence. It is assumed 100 
% detection efficiency, distance to the supernova 10 kpc and target mass 1 kton of ^°^Pb. Upper row: 
equal luminosities in IMH (solid), NMH with small ^13 (dashed). Lower row: L^,^ = 2L^^ (Lp^ = L^^) 
in IMH (solid), NMH with large ^13 (dashed), NMH with small 6*13 (dash-dotted). Left panel: average 
energies from top to bottom {E^J = 13,15 and 18 MeV. Right panel: {E°J — 21 (higher) and 
{E^J = 25 MeV (lower curves). 



Without such an information, HALO alone can give an indication on the presence of energy 
equipartition or on L^^ being larger than L^^ . 

It is clear that to be able to extract the most from future observations, a precise mea- 
surement of neutrino-lead cross sections is called for. This could be realized either at a low 
energy beta-beam facility proposed by [71], or nearby one of the future intense Spallation 
Sources [72, 73]. Indeed a broad physics program can be realized if a low energy neutrino 
facility covering from a few MeV to 100 MeV energy range becomes available in the future, 
including cross section measurements, fundamental interactions tests (neutrino- nucleus co- 
herent scattering, the measurement of the Weinberg angle at low momentum transfer, a CVC 
test, neutrino oscillations to sterile) and measurements of interest for core-collapse supernova 
physics (see [42, 71, 72] and references therein). 

We have summarized our results in figure 11, where all the In-and 2n-event rates are 
plotted in different physical scenarios with a fixed non-electron-type primary neutrino average 
energy. In these figures a given measurement would correspond to just one point. It can be 
seen that if the primary {E^J and {E^^) differ by little, the curves have strong overlap, 
making more difficult to see any evidence of flavor conversion. Obviously, more the primary 
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Figure 9. Relative statistical error ^'^^^"j^^"'^ for the ratio iVi„/A^2n shown in Figure 8 depending 
on the numbers of events. Left: The whole range of possible event numbers in our considerations. 
Right: Zoomed from the left figure showing only for the event numbers < 50. Above each contour 
line the statistical error is smaller than the value indicated in the figure. For instance, assuming 30 
one-neutron and 12 two-neutron events, their ratio is 2.5 and the statistical error about 35% of the 
ratio: ^i„/A^2n = 2.5 ± 0.75. 



average energies differ, easier it is to identify a given physical scenario. In the extreme case 
where the primary {E^^) and {E^J) differ as much as 15 MeV most of the scenarios are very 
likely to be distinguished. 

The whole set of our results for the one- and two-neutron number of events associated 
with a typical cooling phase are shown in figure 12. The straight lines correspond to an 
example where 120 one- neutron events and 30 two- neutron events are measured during the 
explosion, with the associated statistical error. In this example case, one can see that, while 
from the point of view of neutrino properties (mass hierarchy and third neutrino mixing 
angle) all scenarios are possible, rather tight constraints on the primary flux parameters 
- average energy and pinching parameter - can be obtained. Notice that, if the In- and 
2n-event numbers are 210 and 140, respectively (corresponding to {E^J) ~ 25 MeV and 
~ 2 — 4), one can obtain also a clear indication on the mass hierarchy, the value of ^13 
and the luminosity case : the most favorable is L^^ = 2Ly^ and NMH with small ^13). It can 
be also seen that the indication on the presence or absence of the energy equipartition can 
be gained with smaller {E^^) if other information - mass hierarchy and/or ^13 - is known 
from other experiments, as mentioned above. 

In the present work we propagate neutrino probabilities instead of amplitudes. Thus, 
phase effects coming from the shock wave and turbulence effects are not included. Neverthe- 
less, our results represent a good reference of what can be measured. First of all, they give 
an adequate description at early (late) times, when the shock wave has not yet reached (left) 
the MSW region, that will be well identified given the good time resolution of the HALO 



- 22 - 




Figure 10. Ratios of two-neutron event rates as a function of primary non-electron-type neutrino 
pinching parameter a,^^ with different primary non-electron-type neutrino average energies (as in 
figure 6). Upper figures: results of equal luminosities to L^^ ~ 2L^_. in IMH (left) and NMH with 
small (right) . Lower figures: results in IMH to NMH with small 6*13 in the case of equal luminosities 
(left) and L^^ = 2L^^ (right). 

detector. Shock waves and turbulence only impact the signal in the case of normal hierarchy 
and large third neutrino mixing angle when the shock wave has reached the MSW region. 
For such a case one can consider that the neutrino probabilities average approximately to 
1/2, giving a result that is intermediate between the adiabatic and non-adiabatic cases we 
have been calculating. This is also a good approximation for turbulence effects as far as 
the matter density fluctuations keep small enough since large matter density fluctuations 
might break the HL factorization employed in the present work [56, 57]. Obviously, to have 
a detailed description of neutrino flavor conversion and an accurate prediction of the event 
rates in all cases, the inclusion of shock waves and turbulence should be implemented in a 
future calculation. 

4 Conclusions 

We have presented new predictions of the expected neutrino events from iron core-collapse 
supernovae in a lead-based observatory. Such a detector. Helium and Lead Observatory 
(HALO), is currently under construction at SNOLAB. Our predictions are based upon the 
formalism of factorized dynamics. In particular our calculations include collective flavor 
conversion and the MSW effects, as well as decoherence coming from the propagation up to 
the Earth. Determining the luminosities, primary average energies and pinching parameters 
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Figure 11. (Color online) One- and two-neutron emission event rates with different values of the 
primary non-electron neutrino pinching parameter a^^: at the top of each curve a,^^ = 2, at the 
bottom a,j^ = 7. Results include v — ly interactions, the MSW effect and decoherence. The primary 
non-electron neutrino average energies are fixed: (£'" ) = f3 MeV (top left), {E^ ) — 18 MeV (top 
right) and {E^ ) = 25 MeV (bottom figure). Solid lines are for equal luminosities (thick IMH, 
thin NMH with small 613), others for L,,^ = 2L^^: dotted IMH, dashed NMH with large 6*13 and 
dash-dotted NMH with small ^13. The total luminosity is taken to be 10^^ erg s~^, distance to the 
supernova 10 kpc, a target mass 1 kton of ^°*Pb and 100 % detection efficiency is assumed. Notice 
that the results in the case of equal luminosities in NMH with large ^13 are the same as the ones in 
IMH. 



allows an identification of the primary neutrino fluxes at the neutrinosphere, which furnishes 
valuable information to test our understanding of the neutrino transport in the supernova. 

We have shown that in case neutrino flavor transformations occur in supernovae the 
measurement of In- and 2n-event rates as well as of their ratio is particularly sensitive to the 
non-electron neutrino primary average energy and pinching parameter and that interesting 
information can be extracted. Indeed, using information from other detectors, the combi- 
nation of In- and 2n-event rates should allow to identify degenerate solutions of average 
energies and pinching values. Moreover, from the ratio of these events, HALO alone can 
be used to give constraints on these parameters. If the non-electron type primary average 
energy is identified using other detectors one can get even better constraints on the primary 
pinching parameter. 

The actual equipartition of the gravitational energy among the neutrino flavors is also an 
unknown. If the muon or tau neutrino luminosities are larger than the electron neutrino ones 
- as some supernova simulations of the cooling phase seem to indicate - a high energy split 
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Figure 12. (Color online) One- and two-neutron emission event rates with different values of the 
primary non-electron neutrino pinching parameter a^^: at the top of each curve a^^ = 2, at the 
bottom a^^ = 7. Results include i^ — iy interactions, the MSW effect and decoherence. Different colors 
correspond to different non-electron type primary average energies (as in figure 6), solid lines are for 
equal luminosities (thick IMH, thin NMH with small O13), others for L^^ — 2Lj,^: dotted IMH, dashed 
NMH with large 6113 and dash-dotted NMH with small 613. The total time-integrated luminosity is 
taken to be 10^^ erg (~ total events during the cooling phase), distance to the supernova 10 kpc, a 
target mass 1 kton of ^"^Pb and 100 % detection efficiency is assumed. Notice that the results in 
the case of equal luminosities in NMH with large 613 are the same as the ones in IMH. In addition 
is shown an example measurement case with 120 one- neutron and 30 two- neutron events (solid black 
lines) together with estimated statistical errors (dashed black lines). 



might appear, producing different event rates compared to the case where the assumption 
of equipartition of energy among the neutrinos is made. From the combination of In- and 
2n-event rate measurement one can give an indication on the presence or absence of this 
equipartition depending on the actual values of the flux parameters. More the primary 
neutrino average energies differ and if other information (mass hierarchy and/or ^13) is known, 
better indications can be given. 

The present work emphasizes the interest of having more information on the character- 
istics of the high-energy component of the neutrino distributions at the neutrinosphere from 
future supernova simulations: the flux tail is essential to precisely predict the events in such 
a detector. Examples of this type are not only HALO that has been the object of the present 
study but also scintillator detectors using carbon or a Cherenkov detector having oxygen. 
All these have the interesting characteristics of being sensitive to specific features of the high 
energy tail of the neutrino fluxes in a signal from a future supernova explosion. 



- 25 - 



Our specific study for a lead-based observatory furnishes a good example of how, having 
a network of detectors with different energy thresholds, constitute a unique tool to probe 
different components of the neutrino fluxes from a core-collapse supernova explosion and to 
unravel interesting information on the neutrino emission and on neutrino properties. The 
capability of extracting as much information as possible from a future supernova neutrino 
signal in HALO will obviously rely upon the precision with which measurements will be 
made. This involves, in particular, a precise knowledge of the associated cross sections that 
could be at reach with a low energy beta-beam facility or nearby future Spallation Sources. 
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